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We derive the second-order hydrodynamic equation and the microscopic formulae of the relaxation times 
as well as the transport coefficients systematically from the relativistic Boltzmann equation. Our derivation is 
based on a novel development of the renormalization-group method, a powerful reduction theory of dynamical 
systems, which has been applied successfully to derive the non-relativistic second-order hydrodynamic equa¬ 
tion. Our theory nicely gives a compact expression of the deviation of the distribution function in terms of 
the linearized collision operator, which is different from those used as an ansatz in the conventional fourteen- 
moment method. It is confirmed that the resultant microscopic expressions of the transport coefficients coincide 
with those derived in the Chapman-Enskog expansion method. Furthermore, we show that the microscopic 
expressions of the relaxation times have natural and physically plausible forms. We prove that the propagating 
velocities of the fluctuations of the hydrodynamical variables do not exceed the light velocity, and hence our 
second-order equation ensures the desired causality. It is also confirmed that the equilibrium state is stable for 
any perturbation described by our equation. 
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I. INTRODUCTION 

The experiments of relativistic heavy ion collision at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na¬ 
tional Laboratory and the Large Hadron Collider (LHC) at 
CERN seem to have created a hot matter that is most likely 
to be composed of quarks and gluons, i.e., the quark-gluon 
plasma (QGP) Qll. One of the most surprising findings is 
that the created matter is well described by the hydrodynamics 
with tiny dissipation innn. Such a finding prompted an inter¬ 
est in the origin of the viscosity in the gauge theories and also 
the dissipative hydrodynamic equation. The relativistic dis¬ 
sipative hydrodynamic equation is also utilized in analyzing 
various high-energy astrophysical phenomena ifT^ including 
the accelerated expansion of the universe by bulk viscosity of 
dark matter and/or dark energy ifl^ ITsIl . 

One must say, however, that the theory of relativistic hydro¬ 
dynamics for a viscous fluid has not been established on a firm 
ground yet, although there have been many important studies 
since Eckart’s pioneering work ifT^ : A naive relativistic ex¬ 
tension of the Navier-Stokes equation has fundamental prob¬ 
lems such as ambiguity of flow velocity existence 

of unphysical instabilities [|23|, and lack of causality j^l^ . 
the last of which motivated people to introduce the second- 
order hydrodynamic equation. The way of formulation of the 
second-order hydrodynamics is, however, controversial and 
an established equation has not been obtained although some 
suggestive and promising approaches have been proposed. 

It is worth emphasizing here that the second-order hydrody¬ 
namic equation that is free from the causality problem even in 
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the non-relativistic regime is not yet established, either. The 
causality problem inherent in the first-order equation, which 
we call generically the Navier-Stokes equation, appears as the 
instantaneous propagation of information, which is attributed 
to parabolicity of the equation lE^ - l^ . In the seminal pa¬ 
per by Grad he showed that the causality problem could 
be circumvented by the moment method, which now bears 
his name; It is found that the thirteen-moment approximation 
to the functional forms of the distribution function leads to 
a hyperbolic non-relativistic equation satisfying the causality, 
with finite propagation speeds of physical quantities. Here, 
we make a sideremark that the description by the Grad equa¬ 
tion may be called mesoscopic 0^ 1^^ since it occupies an 
intermediate level between the descriptions by the Navier- 
Stokes equation and the Boltzmann equation; see also ll^ . 
It should be noted, however, that the dynamics described by 
the Grad equation has been recently shown inconsistent with 
the underlying Boltzmann equation in the mesoscopic scales 
of space and time ifs^l . Indeed, Grad’s moment method lacks 
a principle for determining the functional form of the distribu¬ 
tion function that is consistent with the underlying Boltzmann 
equation, and then it is inevitable for the moment method to 
adopt an ad-hoc but seemingly plausible ansatz for it. Al¬ 
though there are subsequent attempts to construct the equation 
that respects both of the causality and the consistency with 
the Boltzmann equation in the mesoscopic regime 
the consistency between the resultant equations and the meso¬ 
scopic dynamics of the Boltzmann equation remains unclear. 

Nevertheless, many attempts were made to extend Grad’s 
moment method to establish a mesoscopic description of a 
relativistic system ll^ 1^ [39l - l4^ . but with only a partial 
success, as anticipated. Eor instance, the celebrated Israel- 
Stewart equation lf25l] . which is a typical second-order rela¬ 
tivistic hydrodynamic equation derived from the Boltzmann 
equation based on the moment method with fourteen moments 
employed is found to be incomplete if not incorrect because 
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the solutions behave differently from those of the relativis¬ 
tic Boltzmann equation quantitatively ll39l - |4d]] . The incom¬ 
pleteness or incorrectness can be traced back to ambiguous 
heuristic assumptions inherent in the moment method. Quite 
recently, however, some heuristic but promising methods have 
been proposed 14^ ItM to get rid of such drawbacks, and 
it seems that the resultant solutions indeed become closer to 
that of the Boltzmann equation. Although their results are en¬ 
couraging, one must say that their derivation are still based 
on plausible but ambiguous assumptions that require a mi¬ 
croscopic foundation. In fact, the constitutive equations con¬ 
tain the second-order spatial derivatives of the hydrodynami- 
cal variables, which necessarily leads to the parabolicity that 
should have been avoided. 

Recently, the mesoscopic dynamics or the second-order dy¬ 
namics that respects the causality has been extracted from the 
Boltzmann equation for the non-relativistic case in the clas¬ 
sical regime without recourse to any ansatz for the functional 
forms of the distribution function by two of the present authors 
ll^ : There “the renormalization-group (RG) method’ ’ & 
1^ was adopted as a powerful method of the reduction theory 
of dynamical systems to reduce the Boltzmann equation to the 
mesoscopic dynamics. The basic observation in these works 
is that the first-order hydrodynamics is the slow dynamics 
achieved asymptotically in the kinetic equation 116^ [6411 : The 
asymptotic dynamics is described by the zero modes of the 
linearized collision operator, which happen to be temperature, 
density, and fluid velocity, i.e., the hydrodynamic variables. 
In terms of the reduction theory of dynamical systems ll^ . it 
means that the hydrodynamic variables constitute the natural 
coordinates of the invariant/attractive manifold of the space of 
the distribution function in which the asymptotic dynamics in 
the hydrodynamical regime is described. In the RG method, 
the hydrodynamic variables, which is now the would-be zero 
modes, acquire the time dependence by the RG equation. The 
resultant evolution equation is nothing but the hydrodynamic 
equation, the Navier-Stokes equation. Then, the extension to 
the extraction of the mesoscopic dynamics consists of devel¬ 
oping the way to include some excited (fast) modes properly 
as additional components of the invariant/attractive manifold; 
note that the mesoscopic dynamics is faster than that described 
by the Navier-Stokes equation. In the following natural 
conditions are found to give the adequate excited modes to be 
incorporated in the hydrodynamical variables in the classical 
and non-relativistic case; (A) the resultant dynamics should 
be consistent with the reduced dynamics obtained by employ¬ 
ing only the zero modes in the asymptotic regime; (B) the 
resultant dynamics should be as simple as possible because 
we are interested to reduce the dynamics to a simpler one; the 
term “simple” means that the resultant dynamics is described 
with a fewer number of dynamical variables and is given by 
an equation composed of a fewer number of terms. Here, we 
note that the latter principle (B) is one of the fundamental prin¬ 
ciple of the reduction theory of the dynamics as emphasized 
by Kuramoto if^ . It was shown that these conditions lead 
to a concise scheme called the doublet scheme, and that the 
resultant equation with thirteen dynamical variables satisfies 
the causality in an apparent way and has the same form as that 


of the Grad equation but with different microscopic formulae 
of the transport coefficients and relaxation times; the expres¬ 
sions of the transport coefficients coincide with those by the 
Chapman-Enskog method the novel formulae of the re¬ 
laxation times allow a natural physical interpretation as the 
relaxation times. This is an encouraging result! 

A comment is in order here on the relation between this 
work and if^ in which the two of the present authors (K.T. 
and T.K.) derived a second-order hydrodynamic equation 
from the relativistic Boltzmann equation on the basis of the 
RG method. The derivation presented in Jt^I, however, con¬ 
tained an inconclusive part which is, in retrospect, incor¬ 
rect, unfortunately. Indeed the functional form of the ex¬ 
cited modes was not determined so as to solve the Boltz¬ 
mann equation but that adopted in the Israel-Stewart fourteen- 
moment method was mistakenly used as a possible solution; 
It is known that the Israel-Stewart ansatz does not solve the 
relativistic Boltzmann equation. In this work, we first And a 
proper solution to the relativistic Boltzmann equation in the 
relevant kinetic regime on the basis of an elaborated doublet 
scheme in the RG method, and thus derive the functional form 
of the excited modes that is consistent with the underlying 
Boltzmann equation. Then simply applying the RG equation, 
we obtain the second-order relativistic hydrodynamic equa¬ 
tion, which accordingly gives the correct asymptotic dynam¬ 
ics of the Boltzmann equation in the mesoscopic regime. 

The present paper is an extension of the previous work 14^ 
to a relativistic case with the full quantum statistics as well 
as classical one. We here remark that preliminary results in 
the classical case were announced in IITlI] . Needless to say, 
the quantum statistics is essential in investigating the behav¬ 
ior of a quantum fluid composed of bosons and/or fermions. 
In the present paper, we shall give a detailed and complete 
account of the derivations of the causal hydrodynamic equa¬ 
tions within the quantum and classical statistics together with 
those of the microscopic expressions of the transport coeffi¬ 
cients and relaxation times. We shall also show that a concise 
and natural derivation is possible for the excited modes that 
is given by the doublet scheme on the basis of the very 
principle of the reduction theory of the dynamics. 

Moreover, we prove that the propagating velocities of the 
fluctuations of the hydrodynamical variables do not exceed the 
light velocity, and hence our seconder-order equation ensures 
the causality as desired. It is also shown that the equilibrium 
state is stable for any perturbation described by our equation. 
We give a compact expression of the deviation of the distribu¬ 
tion function to be used in the fourteen-moment method. 

This paper is organized as follows; In Sec. [Ill we briefly 
summarize the basics of the relativistic Boltzmann equation. 
In Sec. imi we derive the causal hydrodynamic equation by 
applying the doublet scheme in the RG method, and give the 
microscopic representations of the transport coefficients and 
relaxation times. Then the basic properties of the resultant 
equation including the causality are shown together with a 
comparison of the microscopic expressions with those given 
by other methods. We devote Sec. |IV]to a summary and con¬ 
cluding remarks. In Appendix we derive the functional 
forms of the excited modes which are given by the faithful so- 
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lution of the Boltzmann equation. In AnnendixlBl the explicit 
solution is given for a linear equation with a time-dependent 
inhomogeneous term appearing in the text. In Appendix O 
we present a detailed and lengthy derivation of the relaxation 
equations, which shows how the microscopic expressions of 
the relaxation times and lengths are obtained. In Appendix iDl 
we give a proof that our second-order equation is really causal 
and that the static solution is stable against any fluctuations. 

In this paper, we use the natural unit, i.e., = c = fee = 1, 

and the Minkowski metric = diag(-|-l, —1, —1, —1). 


metry properties, 

J dp(ppC[f]p = j dpdpidp2dp3w(p, pi|p2 , Ps) 

X (<Pp -f Ppi — Pp2 ~ fPpa) 

X ((1 + a/p)(l -f afp,)fp2fp2 
-/p/pi(l + a/p2)(l + a/p3))- (6) 

Substituting (l,p^) into (pp in Eq. (|6|, we find that are 

collision invariants satisfying 


II. PRELIMINARY 

In this section, we summarize the basic facts about the rel¬ 
ativistic Boltzmann equation 


A. Basics of relativistic Boltzmann equation 

The relativistic Boltzmann equation reads if^ 1?^ 

P^dpfpix) = C[f]p{x), (1) 

where fp{x) denotes the one-particle distribution function 
with being the four-momentum of the on-shell particle, 
i.e., p^Pp = p^ = m? and p° > 0. The right-hand-side 
term C[f]p{x) denotes the collision integral 


/dp^C[/]p = 0, (7) 

J = 0. (8) 

due to the particle-number and energy-momentum conserva¬ 
tion in the collision process, respectively. We note that the 
function ipop = Q!(a;) -|-p^/3p(a;) is also a collision invariant 
where a{x) and /3^(x) are arbitrary functions of x. 

Owing to the particle-number and energy-momentum con¬ 
servation in the collision process leading to Eqs. (|7]i and ®, 
we have the balance equations 

= 0, (9) 

where the particle current and the energy-momentum ten¬ 
sor are defined by 


C{f]p{x) - dpidp2dp3a;(p,pi|p2,p3) 

m = J dppVp, 

(11) 

X + afp{x)){l + afp,{x))fp2{x)fp2{x) 

- fpix)fp,ix){l + afp2ix)){l + afp^ix))), (2) 

= J dpp^Vp, 

(12) 


where w(p,pi|p 2 ,p 3 ) is the transition probability due to the 
microscopic two-particle interaction with the symmetry prop¬ 
erty 

W(P,P1|P2,P3) =0o{p2,P3\P,Pl) 

= w(pi,p|p3,P2) = a;(p3,p2|pi,p), (3) 

and the energy-momentum conservation 


respectively. It should be noted that any dynamical properties 
are not contained in these equations unless the evolution of fp 
has been obtained as a solution to Eq. O- 

In the Boltzmann theory, the entropy current may be de¬ 
fined iF^ by 


/ rlnnM 

J- i„ V- (1 + “/p) Ml + afp) 

/ ^PP 



w(p,Pl|p2,P3) oc S'^ip+Pi -P 2 -Pa), 


(4) 


The entropy current satisfies the divergence equation 


and a represents the quantum statistical effect, i.e., a = +1 for 
boson, a = — 1 for fermion, and a = 0 for the Boltzmann gas. 
In the following, we suppress the arguments x, and abbreviate 
an integration measure as 

dp EE d^p/[(27r)V], (5) 

with p being the spatial components of the four momentum 
p^ when no misunderstanding is expected. 

Eor an arbitrary vector pp iH, the collision integral satis¬ 
fies the following identity thanks to the above-mentioned sym- 


= - 


dpC[f]p\n 


fp 

1 + 0.fp 


(14) 


because of Eq. O- One sees that is conserved only 
if ln(/p/(l + afp)) is a collision invariant, i.e., ln(/p/(l + 
afS ) = Pop = a{x) + p^l3p{x). One thus finds iIt^ 
1^ that the entropy-conserving distribution function can be 
parametrized as 


fp — 


1 


feq 

Jv ’ 


(15) 
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where T, fi, and may depend on the space and time, and 
are interpreted as the local temperature, chemical potential, 
and flow velocity, respectively, with the normalization 

u^^Ufj, = 1. (16) 

Thus the function ( fTSl l is identified with the local equilibrium 
distribution function. We see that the collision integral identi¬ 
cally vanishes for the local equilibrium distribution /p'^ as 

cn = 0 , ( 17 ) 

owning to the detailed balance 

io{p,Pi\P2,P3) [(1 + af;^)ii + 

- /;Vp7(i + «/p7)(i + «/p7)] = (18) 

guaranteed by the energy-momentum conservation (01. 
Substituting fp = into Eqs. (fTTT i and (fTSl i. we have 

= nu^^ = N^, (19) 

= eu^^u'' - PA'"'" = Tg'"'", (20) 

with 

A'"'" = 5'"'"(21) 

Here, n, e, and P denote the particle-number density, internal 
energy, and pressure, respectively, whose microscopic repre¬ 
sentations are given by 


n = y dpf^^{p ■ u) 

OO 

= (27r)“^47rm^ {k'm/T)~^K 2 {km/T), 

fe=i 

( 22 ) 

e = j dp/ 7 (p ■ uf 

[km/T)~^K^{km/T) 

— 77777 _ — _ - _ 

EZi a'-ie''"/^(Zm/T)-iiT2((m/T) 

EZ=1 a’^-^e>^>^/^{km/T)-^K2{km/T)'\ 

P = J 

with K 2 {z) and K'i{z) being the second- and third-order mod¬ 
ified Bessel functions. Setting a = 0 in the above expressions, 
we can check that the classical expressions for n, e, and P if^ 
are reproduced. We note that and in Eqs. (fTOl l and 
(l20l l are identical to those in the relativistic Euler equation, 
which describes the fluid dynamics without dissipative effects, 
and n, e, and P defined by Eqs. (l22]) - (l24l) are the equations 
of state of the dilute gas. Since the entropy-conserving distri¬ 
bution function reproduces the relativistic Euler equation, 
we find that the dissipative effects are attributable to the devi¬ 
ation of fp from fp‘^. 


III. RELATIVISTIC CAUSAL HYDRODYNAMICS BY 
DOUBLET SCHEME IN RG METHOD 


In this section, we derive the causal relativistic hydrody¬ 
namic equation as the mesoscopic dynamics from the rel¬ 
ativistic Boltzmann equation (l27T i: The derivation is based 
on the doublet scheme in the RG method developed for the 
non-relativistic case in ; the present formulation is an ex¬ 
tension to the relativistic case and given in a simplified and 
more transparent manner. We examine some properties of 
the resultant equation, concerning the frame, the stability of 
the equilibrium state, and the causality as well as the micro¬ 
scopic representations of the transport coefficients and the re¬ 
laxation times. It will be noted that our formalism solving 
the Boltzmann equation gives the compact expression of the 
perturbed distribution function, which may be used in the mo¬ 
ment method as the proper ansatz of the distribution function 
that can lead to the hydrodynamic equation consistent with the 
Boltzmann equation. 


A. Reduced dynamics by RG method 

1. Macroscopic-frame vector 

Since we are interested in the hydrodynamic regime to be 
realized asymptotically where the time and space dependence 
of the physical quantities are small, we try to solve Eq. (ID in 
the situation where the space-time variation of fp{x) is small 
and the space-time scales are coarse-grained from those in the 
kinetic regime. To make a coarse graining with the Lorentz 
covariance being retained, we introduce a time-like Lorentz 
vector denoted by o'" with a2 > 0 and aO > 0 i?! [tI, 
which may depend on x'"; o'" = a^{x). Thus, speci¬ 
fies the covariant but macroscopic coordinate system where 
the local rest frame of the flow velocity and/or the flow ve¬ 
locity itself are defined: Since such a coordinate system is 
called frame, we call o'" the macroscopic frame vector. In 
fact, with the use of a^, we define the covariant and macro¬ 
scopic coordinate system (r, cr'") from the space-time coordi¬ 
nate x'" as dr = a'^dx^ and dtr'" = (p'"'" — a^aPjcff)dxv, 
which lead to derivatives given by djdr = {a^ jdf^dp and 
dfdap = (p'"'" — a^aP j df)dy. 

Then, the relativistic Boltzmann equation ([Til in the new 
coordinate system (r, ct'") is written as 

p-a{T,a)-^fp{T,a) +P^^/p(ao') = C'[/]p(T, cr), 

(25) 

where a^{T,a) = a^{x) and /p(r, u) = fp{x). We remark 
the prefactor of the time derivative is a Lorentz scalar and pos¬ 
itive definite; p-a(T, cr) > 0, which is easily verified by taking 
the rest frame of p°. 

Since we are interested in a hydrodynamic solution to Eq. 
(l25l l as mentioned above, we suppose that the time variation 
of a^{T, cr) is much smaller than that of the microscopic pro- 
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cesses and hence a^(r, a) has no r dependence, i.e., lows; 


a^{T,a) = a^{a). (26) 

Then, with the use of Eq. ( l26l l, we shall convert Eq. ( |25] ) into 




' (27) 


~p-a{a) 9(7^ ■ 


/p(r,cr;To) = /®(T,cr;To) + e/^^^(T,cr;ro) 

+ e^/p^^('r,CT;To) H-, (29) 

/p(to, ct) = a) + e4^)(To, cr) 

+ e^/p^^('fo,o-) H-. (30) 

The respective initial conditions at r = tq are set up as 

/^'^(To,cr;To) = 4 ')(to,ct), Z = 0, (31) 


Here, the parameter e is introduced for characterizing the 
smallness of the inhomogeneity of the distribution function, 
which may be identified with the ratio of the mean free path 
over the characteristic macroscopic length, i.e., the Knudsen 
number. Since e appears in front of the second term of the 
right-hand side of Eq. (l27T i. the relativistic Boltzmann equa¬ 
tion has a form to which the perturbative expansion is appli¬ 
cable. 

In the present analysis based on the RG method, the per¬ 
turbative expansion of the distribution function with respect 
to e is first performed with the zeroth-order solution being the 
local equilibrium one, which has no dissipative effects. The 
dissipative effects are taken into account in the higher orders; 
the spatial inhomogeneity as the perturbation gives rise to a 
deformation of the distribution function which is responsible 
for the dissipative effects. Note that the deformation also can 
trigger a relaxation toward the local equilibrium state. Thus, 
the above rewrite of the equation with e reflects the physi¬ 
cal assumption that only the spatial inhomogeneity plays dual 
roles as the origin of the dissipation and the cause of a relax¬ 
ation to the local equilibrium state. It is noteworthy that our 
RG method applied to the non-relativistic Boltzmann equa¬ 
tion with the corresponding assumption successfully leads to 
the non-relativistic causal hydrodynamic equation which 
means that the present approach is simply a relativistic gener¬ 
alization of the non-relativistic case. 


2. Construction of approximate solution around arbitrary initial 

time 

In accordance with the general formulation of the RG 
method if^ 1^ [s^ . let fp{T, a) be an exact solution yet to 
be obtained with an initial condition set up, say at r = —c>o. 
Then we pick up an arbitrary time t = tq in the (asymptotic) 
hydrodynamic regime, and try to obtain the perturbative so¬ 
lution fp to Eq. dZTl l around the time r = tq with the initial 
condition 


In the expansion, the zeroth-order value it;tq) = 

fp^\To,a) is supposed to be as close as possible to an 
exact solution fp{T,a). In the RG method, the globally 
valid solution is constructed by patching the local solutions 
ct; To) which are only valid around r = tq, which is 
tantamount to making an envelope curve of the perturbative 
solutions with tq being t he p arameter characterizing the per¬ 
turbative trajectories 11^ 1^ . 

Substituting the above expansions into Eq. (l271 i we ob¬ 
tain the series of the perturbative equations with respect to 
e, where the macroscopic frame vector is now replaced by a 
T-independent but tq- dependent one 

a'"(CT) = a'"(CT;To). (32) 

We have now a hierarchy of equations in order by order of 
e. As is mentioned before, our strategy to obtain the meso¬ 
scopic dynamics is constructing it as a minimal extension of 
the hydrodynamic one that is to be realized asymptotically af¬ 
ter a long time within the Boltzmann equation. Notice that the 
hydrodynamics is a closed slow dynamics described solely by 
the would-be zero modes of the linearized collision operator 
corresponding to the conservation laws. The slowest dynam¬ 
ics will be given as a stationary solution, which actually exists 
for the zeroth order equation; the stationary solution is noth¬ 
ing but the local equilibrium one ItM It^ . In our way of the 
solution of the Boltzmann equation on the perturbation the¬ 
ory with the single expansion parameter e, the deviation of the 
distribution function from the local equilibrium one is caused 
by the spatial inhomogeneity as given by the perturbative term 
in Eq. dZTl l and hence is proportional to e. We shall show that 
this setting of the analysis successfully solves the Boltzmann 
equation in a consistent way and leads to the mesoscopic dy¬ 
namics. 

With the above order counting in mind, let us construct the 
perturbative solution in the asymptotic regime order by order. 
The zeroth-order equation reads 

■^fp°HT,a;To) = - ^ --C[/(°)]p(r,CT;ro). (33) 

dr ^ p- a{a]To) 


~ \ r f \ Since we are interested in the slow motion which would be 

fp\T — tqjCTjTq) — /p(to,ct), (28) realized asymptotically as r —>■ oo, we should take the fol¬ 

lowing stationary solution, 

where we have made explicit that the solution has the tq de¬ 
pendence. The initial value or the exact solution as well as 
the perturbative solution are expanded with respect to e as fol- 


= 0 , 


(34) 
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which is realized when (r, u; tq) is the fixed point, 
C[f°'>]p{T,a-,To) =0, 


p-a{a;To) 


(35) 


for arbitrary a. We see that Eq. dlST l is identical to Eq. (fTTl) . 
and hence fp^\T, tr; tq) is found to be the local equilibrium 
distribution function (ffSl) : 


fi°\r,a;To) = /p"'(cr;ro) 


(36) 


e[p'"«p(cr;To)-/i(<T;ro)]/T((T;ro) _ q’ 

with u^{a; Tq)u^{(J] tq) = 1, which implies that 

fp°H'ro,cr) = = ro,cr;To) = /p'’(cr;ro). (37) 

The five would-be integral constants T((t;to), /i(cr;To), and 
Mp(cr; To) are independent of t but may depend on tq as well 
as a, and the local temperature, local chemical potential, and 
flow velocity can be naturally obtained. Eor the sake of the 
convenience, we define the following quantity: 


To) = 1 + a/p‘i(cr; To) 


p 

gb''«p(o-;To)-/i(o';To)]/T((T;To) 
e[p'^iip(CT;To)-/j(cr;To)]/T(o-;To) _ ' 


(38) 


We remark about an explicit form of a^{a\TQ) that should 
be a Lorentz four vector described by the hydrodynamic vari¬ 
ables T((t; To), To), and tq) and their derivatives. 
In the case of the first-order hydrodynamic equation, it was 
shown that as long as such a a^ \a] tq) is independent of 
the momentum the leading terms of the resultant equation 
perfectly agree with those obtained with the choice 


a'"{(J-,To) = u'"(cr; To). 


In the present work, we will present the analysis that is based 
on this choice, and derive the second-order hydrodynamic 
equation as a natural extension of the first-order one obtained 
in izl. In the following, we suppress the coordinate argu¬ 
ments (cr; To) when no misunderstanding is expected. 

The choice leads to the following identities 


— = u^d 

dr ^ 

oa,. 


(40) 

(41) 


Note that d/dr and are the Lorentz-covariant temporal 
and spacial derivatives, respectively. 

Now that the preliminary set up is over, let us move to the 
analysis of the first-order equation. Inserting the expansion 
( |29] | into Eq. (l27l i with the setting ( |39] |. we have the first-order 
equation as 

+ /pVr^Op, (42) 


where Lpg is the linearized collision operator 


f = I 

^pq — \Jp Jp ) 


1 1 


1 


p-uSfc 


■C[fh 


f eq req 
J a J a 


2\p ' U J 

£eq £eq £eq 
Jpi Jp2 Jp3 
neq 

Jp 


[ dpidp2dp3Uj{p,pi\p2,P3) 
{^pq + ^Piq ~ Sp2q — Sp^q), 


and Fqp is an inhomogeneous term 


1 




p • u 


(43) 


(44) 


Eor the sake of simplicity, we rewrite Eg. (1421) in a vector form 

= nnknnr^f^^Hr) + nnPo, 

(45) 

where we have treated /®'t and as a diagonal matrix. 

The linearized collision operator has some remarkable 
properties that play important roles in the following analy¬ 
sis. To see this, let us define an inner product for two arbitrary 
functions ipp and Xp by 


(V’,X)= / dp(p-u)/®Vr'*ApXp- 


(46) 


This inner product is a generalization of the one introduced in 
SS for the classical statistics to the quantum one. This 
inner product respects the positive definiteness as 

{ip,ip)>0, ioiipp^O, (47) 

because (p ■ u) in the inner product is positive-definite. Then 
( 39 ) we find that L is self-adjoint with respect to this inner product 


Ilf 

= / dpdpidp2dp3w(p,pi|p2,p3) 

X fPfpkpkpk^P + V'pi - V'P2 - i’ps) 

X iXp + Xpi - Xp2 - Xps) 

= {U,x), (48) 

and non-positive definite 

Ilf 

/ dpdpidp2dp3a;(p,pi|p2,p3) 

X fkfpkpkpk^P + V'pi - V'P2 - V'ps)^ 

< 0, (49) 

with Tpp and Xp being arbitrary vectors. The operator L has 
the five eigenvectors belonging to the zero eigenvalue; 


=0, 


with 


‘Pop = 


a = p, 

1, a = 4. 


(50) 


(51) 
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We note that (pQp with a = 0, ■ ■ ■ , 4 are the collision invari¬ 



ants, and span the kernel of L. We call ipop the zero modes. 

To represent the solution to the first-order equation (l42l) in a 

Qo 

Qi 

comprehensive way, we define the projection operator Pg onto 


the kernel of L which is called the Pg space and the projection 
operator Qg onto the Qg space complement to the Pg space: 


Pi 

[Po V’] p = ^opVo^is (<Po, V') , (52) 

Qo = 1 - Po, (53) 


Po 


where is the inverse matrix of the the P-space metric 
matrix 77 Q ^ defined by 

(54) 

Now the solution to (l42l) is given in terms of Po and Qo as 


/W(T,a;ro)=rV"‘i 




^' + (t - to)PoFo 


+ (ei(— 0 ) _ i)L-^Q^Fo 


(55) 


with 


= PHr = TO,a- to) = (56) 

where is the integral constant. Here, the second and third 
terms in Eq. (fSSl) describe the motion caused by the pertur¬ 
bation term Fq, i.e., the spatial inhomogeneity, while the first 
term can be identified with the deviation from the stationary 
solution /®'t, which should be constructed in the perturba¬ 
tive expansion with respect to the ratio of the deviation from 

/eq {Q /eq Jjj jjjg /eq j g ^ 

is nothing but the time-dependent 
solution to Eq. (l3?t valid up to 0{e). It is obvious that this so¬ 
lution relaxes to + ef^’^Po'i’) in the asymptotic regime, 
because vanishes as T — 00 . In order 

to obtain the time-dependent solution that describes the re¬ 
laxation process to the stationary solution we must sup¬ 
pose that Po4' = 0, i.e., 'P contains no zero modes. This is 
a kind of the matching condition. Indeed, if T* were to con¬ 
tain zero modes, such zero modes could be eliminated by the 
redefinition of the zeroth-order initial value specified by the 
local temperature T((t;to), chemical potential fi{a;To), and 
flow velocity u^{a; tq). In fact, + j5) 

can be written as a sum of the derivatives of /®‘^ with respect 
to T, pL, and Up with the identification, = S{up/T) = 
Sup/T+UpS{l/T) and/3 = -S{fi/T) = -Spi/T-fj.S{l/T), 
which leads to = —fp'^fp{p^^6{up/T) — d{p/T)). We 
note that the transverse component of is proportional to 
Sup. Thus we see that the possible existence of the zero modes 
in T' would be renormalized into the local temperature, chem¬ 
ical potential, and flow velocity, and absorbed into the redefi¬ 
nition of the initial distribution function at local equilibrium. 

We note the appearance of the secular term proportional to 
r—Tg in Eq. ( |55] ). which apparently invalidate the perturbative 
solution when |t — ro| becomes large. 


FIG. 1: Decomposition of the solution space of the Boltzmann equa¬ 
tion. The Po space is the kernel of the linearized collision operator, 
while the Qo space is spanned by excited mode, which is decom¬ 
posed into the Pi and Qi spaces. 


Eor later convenience, let us expand with respect 

to T — Tg and retain the terms up to the first order as. 


P\T,a;To)^nnh + iT 


Tg)LT' + {t - to)PoFo 


+ {t - to)QoFo 


(57) 


Here the neglected terms of 0((r — tq)^) are irrelevant when 
we impose the RG equation, which can be identified with the 
envelope equation 11 ^ and thus the global solution is con¬ 
structed by patching the tangent line of the perturbative solu¬ 
tion at the arbitrary initial time t — To, as mentioned before. 

Now the problem is how to extend the vector space beyond 
that spanned by the zero modes to accommodate the excited 
modes that are responsible for the mesoscopic dynamics and 
should consist of the basic variables together with the zero 
modes to describes the second-order hydrodynamics. The 
vector space to which the excited modes belong are called the 
Pi space. Here one should note that the Pi space is a sub¬ 
space of the Qg space, as shown in Pig. [T] To this end, let 
us see what the first-order solution (l57l i tells us how to extend 
the vector space. In fact, to do that we only have to make the 
following requirement; The tangent spaces of the perturbative 
solution at T = Tg become as small as possible to simplify 
the obtained equation. Instead of the two requirements (A) 
and (B) introduced in Sec. |I] we utilize here this one require¬ 
ment to determine explicit forms of the vector T' and the Pi 
space. We note that although the resultant forms of them are 
the same as those obtained with (A) and (B), the derivation of 
them becomes more natural and straightforward. Simplicity 
of the obtained equation is one of the basic principles in the 
reduction theory of dynamical systems. Here, we note that 
such tangent spaces are spanned by the terms proportional to 
r — To in Eq. ( fSTT i. while the Pi space is spanned by all the 
terms except for the zero modes in Eq. dSTl i. 

Thus, we can reduce this requirement to the following two 
conditions; 

• LT' and Qo^o should belong to a common vector space. 


• The Pi space is spanned by independent components of 
LT' and T'. 
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The first condition is restated as that T' and L should 

belong to a common vector space. Therefore let us calculate 
L~^QqFq and examine the structure of the vector space to 
which it belongs. This explicit calculation of the deforma¬ 
tion of the distribution function constitutes one of the central 
parts of the present work, contrasting to the moment method 
in which some seemingly plausible ansatz is adopted without 
any explicit solution. A straightforward but somewhat tedious 
calculation of it is worked out in Appendix The result is 
given as 


[L ^QoFo]^ - - 






(58) 


Here, Hp, J^, and tt^'' are microscopic representations of dis¬ 
sipative currents whose definitions are given by 


= —(np,jp^7r^‘^), (59) 


with 


'l dP 


3 de 

n 


np = {p-uf ^ ^ +{p-u) 


dP 

dn 




- (60) 

(61) 

(62) 


In the above equations, we have introduced the enthalpy per 
particle h and the projection matrix given by 


h = {e + P)/n, (63) 

^ -P - 2/?>AP''AP‘^), (64) 


Thus we now see that 4/ can be written as a linear combi¬ 
nation of these bases as 




(ff,L-iff) 


n 


h[L-^>]p 




L 3 




J„ 


(69) 


Here we have introduced the following nine vectors as mere 
coefficients of the basis vectors; 


n(cr;ro), JP{a;To), {a;To). (70) 

We stress that the form of T' given in Eq. (1691 ) is the most 
generic expression that makes L'k and QqFq belong to the 
common space. 

As is clear now, we see that the Pi space is identified 
with the vector space spanned by Hp, Jp, [L“^1I]^, 

and The sets of fl and L“^ff, Jp and 

L~^JP, and and L~^7tP'’^ are called the doublet modes 
The Qo space is now decomposed into the Pi space 
spanned by the doublet modes and the Qi space which is the 
complement to the Pq and Pi spaces. The corresponding pro¬ 
jection operators are denoted as Pi and Qi, respectively. 

Now we find that the coefficients JP and in Eq. (l69b 
are taken to be transverse without loss of generality; i.e., 

JP = (71) 

= AP'^P’^TVp^, (72) 

because of Eqs. ( l67l l and (|68]l. The properties ( fTTI ) and (l72b 
lead to the following identities; 

UpJP = UpirP’' = Ap^irP'' = 0, (73) 

7rP‘'=Tr’'P. (74) 


respectively. It is notable that the coefficients of the nine vec¬ 
tors [L“^n]^, [L“^ J^]^, and [i/“^7r^‘"]^ are linearly inde¬ 
pendent, i.e., the following statement is true; 


It will be found that H, Jp, and can be identified with the 
bulk pressure, thermal flux, and stress pressure, respectively. 
The second-order equation is written as 


a(_T-iV • u) + pP{-h-^S/p{p/T)) 

+ r^T-^Ap,p,VPu'^) =0, V T, p, uP 

a = pp = jP’' = 0. (65) 


= nniinny^f^^Hr) 

+ nnK{T-TQ), (75) 


Thus, we can take the following nine vectors 


with the time-dependent inhomogeneous term given by 


[L-iH]^, [L-^J%, ( 66 ) 

as a set of the bases of the vector space that [L^^Qq^o] ^ 
hence T' belong to. Here we note that the above Lorentz vector 
and the tensor are transverse; 

[L-i>]^ = A^‘^[L-ij.]^, (67) 

[L-^nP’']^ = AP-'P^[L-%,]^. ( 68 ) 


k{t - to) = fiPHt) + ^B[{nn)-^p\T))] ' 
= Finn + (r - to)PoFo 


+ (e^(—o) _ i)nQoFo 


2 ^ 


To)^ + {t - to)PoFo 


+ (e 


i(T-To) _ 


l)L-noFo 


(76) 
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Here, Fi and B are matrices and their components are given 
by 


y p P . 

In Eq. (Illll, we have used the notation 


ill) 

req Jeq xeq leq 
J q J q Jr Jr ’ 


/=/■*'» 


(78) 


1. In the present approach, we are solving the Boltzmann 
equation ([T]i as faithfully as possible, in contrast to the 
Israel-Stewart fourteen-moment method 12^ . in which 
an ansatz for the solution is imposed in the form / = 

/eq_^gjeqJeq^l4M with = a-f -f 

Here, the coefficients a, b^, and c^’’' are definite func¬ 
tions of T, /i, u^, n, and JH]. It is interesting 

that our initial value given in Eq. ( |69] | provides a 
foundation of the fourteen-moment method but with a 
novel form of vj/i^M 


[BiI^'x\^ = j AqArBpqrlijqXv (79) 

The solution to Eq. ( fTSl ) around t ^ tq is found to take the 
following form 

p\T,a-To) = nn[{T-To)p^ 

+ {t - to){L - d/ds)Pig{s)Qo 

- (1 + (t - To)d/ds)Qig{s)Qo\K{s) , 

J s=0 

(80) 


the initial value of which reads 


f ^^\ a ; To ) = = to,ct;to) 

= -nnQig{s)QoK{s) 


( 81 ) 


The derivation of this solution is presented in Appendix |Bl 
where the complete expression of the solution not restricted 
to r ~ To is given: In Eq. dSOl l. we have retained only terms 
up to the first order of (r — tq), and introduced a “propagator” 
defined by 


g{s) = {L-d/ds)-\ (82) 


We notice again the appearance of secular terms in Eq. dSOl l. 

Summing up the perturbative solutions up to the second or¬ 
der with respect to e, we have the full expression of the ap¬ 
proximate solution around r ~ tq to the second order: 


/(t, cr; To) 

= n + enn [(1 + (r - ro)m + {t- to)Fo 
+ - To)Po + (r - ro)(L - d/ds)Pig{s)Qo 


- {1 +{t - To)d/ds)Qig{s)Qo K{s) 


s=0 


(83) 


with the initial value 


f{a-,To) = n + enn'if 

- e^nnQig{s)QoK{s) 


(84) 


We note that the possible appearance of the fast motion caused 
by the Qi space in Eq. dS^ is avoided by an appropriate 
choice of the initial value ( l84b . as in the first-order solution; 
see AppendixlBlfor the detail. 

A couple of remarks are in order here: 


2. Expanding g{s)Qo in terms of L~^d/ds, the term 
(/(s)Qo7f(s)|s=o in Eqs. ( |83] ) and ( l84b is reduced to 
the form of infinite series as 


g{s)QoK{s) 


n —0 


(85) 


because d"'K{s)/ds^\s=o does not vanish for any n\ 
see Eq. ( |76] |. Admittedly the existence of such an infi¬ 
nite number of terms would be undesirable for the con¬ 
struction of the (closed) mesoscopic dynamics. It will 
be found, however, that an averaging procedure for ob¬ 
taining the mesoscopic dynamics nicely leads to a can¬ 
cellation of all the terms but single term in the resultant 
equation of motion thanks to the self-adjointness of L 
and the structure of the Pi space spanned by the doublet 
modes; see Eq. ( l94b below. 


3. RG improvement of perturbative expansion 


We note that the solution dS^ contains secular terms that 
apparently invalidate the perturbative expansion for t away 
from the initial time tq. The point of the RG method lies 
in the fact that we can utilize the secular terms to obtain an 
asymptotic solution valid in a global domain. Now we see that 
/p(r, cr; To) in Eq. d83]) provides a family of curves parameter¬ 
ized with Tq. They are all on the exact solution fp{a; r) given 
by Eq. d84l) at r = tq up to O(e^), but only valid locally for r 
near tq. Thus, it is conceivable that the envelope of the fam¬ 
ily of curves, which is in contact with each local solution at 
r = Tq, will give a global solution in our asymptotic situation 
0^ 1^ 1^ 1^ |6J]. According to the classical theory of en¬ 
velopes, the envelope that is in contact with any curve in the 
family at t = tq is obtained by 




= 0 , 


( 86 ) 


where the subscript p is restored for later convenience. Equa- 
tion d86] ) is called the renormalization group equation 114^ . 
and has also the meaning of the envelope equation ifsill . We 
call Ea.d86ll the RG/Envelope or RG/E equation following 
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. Now Eg. (1861 1 is really reduced to 
dr 


.^i + en^))-enn 


L'$ + PqFq + QqFq 




Po + {L- d/ds)Pigis)Qo 


-{d/ds)Qig{s)Qo 


K{s) +O{e^)=0. (87) 

s=0 


It is noted that Eq. dSTl i gives the equation of motion govern¬ 
ing the dynamics of the would-be fourteen integral constants 
T(tT; r), ^(ct; r), u^{a] r), n(CT; r), r), and r). 

The envelope function is given by the initial value (l84l l with 
the replacement of tq = r as 


= fp{T,a]To = t) 
= /p(ct;to = r) 

= mi + erm 


- e^mmiGis)QoK{s) 


+ 0{e^), 

ro=T 

( 88 ) 


where the exact solution to the RG/E equation (IStI i is to be in¬ 
serted. We note that the envelope function (t, a) is actually 
the global solution that solves the Boltzmann equation dZTl i up 
to 0{e^) in a global domain in the asymptotic regime: Indeed, 
for arbitrary r(= tq) in the global domain in the asymptotic 
regime, we have 


^fp (r,a) = —fp(T,a;To) 

d f , . 

= ^/p(go-;to) 


To=r 


d ~ 

+ ■^fp{T,cr-,To) 
OTq 


To=T 

(89) 


where the RG/E equation dShl) has been used. Eurthermore, 
since /p(r, u; tq) solves Eq. dZTl i with a^(cr) = M^(cr; tq) up 
to O(e^), the r.h.s. of Eq. (18^ reads 


d ^ ^ _^ _ 1 

OT p-u{a;To) 


1 


— e 


p-u{a;To)^ 


C[f]p{'r,cr;To) 

(90) 


Then inserting the definition of f^{T, a) given in the first line 
of Eq. dSST l. we have 


-/f (tct) = 


1 


:C'[/°]p(go-) 


p - u { a - T ) 

-e-7--f (9(6^). (91) 

p • u{a; T) 0(7^ ^ 


This concludes the proof that the envelope function (r, a) 
is the global solution to the Boltzmann equation dZTl i up to 
O(e^) in a global domain. 


It is noteworthy that we have derived the mesoscopic dy¬ 
namics of the relativistic Boltzmann equation ( l27b in the form 
of the pair of Eqs. (|87] | and (ISSI l. It is to be noted that an infi¬ 
nite number of terms, produced by g{s), are included both in 
the RG/E equation and the envelope function. 

We observe that the RG/E equation ( l87T i includes fast 
modes that should not be identified as the hydrodynamic 
modes even in the second order ones. While these modes 
could be incorporated to make a Langevnized hydrodynamic 
equation, we average out them to have the genuine hydrody¬ 
namic equation in the second order. This averaging can be 
made by taking the inner product of Eq. dSTl l with the zero 
modes vJqj, and the excited modes [L“^(II, used 

in the definition of T'p. The first averaging leads to 


f 

', , d 


j ‘iPV’op 

(j> ■ u)— +ep-V 

fP{l + e/pm 


= 0 + O(e3), (92) 

and the second averaging 


J dp[L-\fi,j>^m)]j 





= e{L-Hfl, >,77^*^), XT') 

+ + 0{€^). 


Here we have used the identity given by 


(93) 


(T-i(fl, >, (L - d/ds)Pig{s)QoKis) ) 

s—0 

= {{L - d/ds)L-\U, m,PiG{s)QoK{s) ) 

= {{L - d/ds)L-\U, m,Gis)QoK{s) ) 

s=0 

= {L-\ll,r,m,iL-d/ds)gis)QoK{s) ) 

s=0 

s=0 

= (L-i(n,>,7f'^"),iT(0)) 

+ (94) 


where utilized are the self-adjointness of Lpq shown in Eq. 
(l48l l. the structure of the Pi space spanned by the dou¬ 
blet modes, i.e., the pairs of lip and [L“^Il]^, and 

[L“^ and TT^'^ and and the equality iT(0) = 

J 7 i/eq /2 derived from Eq. (|76ll. 

Thus the pair of Eqs. ( |92] | and (|93] | constitutes the hydro- 
dynamic equation in the second order, i.e., the equation of 
motion governing T, /i, u^, 11, J^, and ir^''. It is to be noted 
that this pair of equations is free from an infinite number of 





























11 


terms in contrast to the RG/E equation dSTl i and much simpler 
than it. We stress that this simplification through the averag¬ 
ing by L“^(n, J^, TT^^) is due to the self-adjointness of L and 
the structure of the Pi space spanned by the doublet modes 
andL-i(n,>,7r'''^). 


B. Properties of the reduced dynamics 


We now put back to e = 1. Noting that {p ■ u)-^ p -V = 
p^dfj,, we find that Eq. (l92l i finally takes the following form 


« Ttia _ Q 


(95) 


with 


J^^a _ 

hydro — 


- (P + a = v, 


nu^ 




a = 4. 


(96) 


We remark that Eq. (I^l l is nothing but the balance equa¬ 
tions and and can be identified with the energy- 

momentum tensor and particle current iV^ in the Landau- 

Lifshitz frame, respectively. Indeed, we can derive the same 
expression as by substituting the distribution function 

/^(r, a) in Eq. dSST l into the definitions of and given 
by Eqs. (fT^ and (fTTT l. 

After a straightforward manipulation whose details are pre¬ 
sented in AppendixlO we can reduce Eq. (|9^ into the follow¬ 
ing relaxation equations: 


n = -C 6 » 

Q 

— Tn^n — £njV • J 

OT 

+ Knnll^ 

+ fennll^ + bnjjJ^’Jp 
u 


^nTTTr”^^ 




- TjA^^P^Jp - f jnV^n - 
OT 


(97) 


+ K^\^nv^r ■ 




+ + KfjJpa^^P + nfjJpUj^^P 


''JJ'- 

.( 1 ), 


+ bjYij^J^ + bjjT^JpTT^^, 


A* 


T 


(98) 




+ + b^jjj‘^^r'> + (99) 


where we have introduced the notation Ap^ 

for a traceless and symmetric tensor. Here 6 = V • m, 
= AP’^P'^Vp^, and - V^u'') denote 

the scalar expansion, shear tensor and vorticity, respectively. 
We refer to Appendix |C] for the explicit definitions of many 
other hydrodynamic valuables introduced in (l97li-(l99ll. 

Now the physical meaning of each term in (l97]) - (l9^ should 
be clear: The first lines in Eqs. (|97]|-(|99]| are identical with 
the so-called constitutive equations, which define the rela¬ 
tions between the dissipative variables H, Jp, and ttP'' and 
the thermodynamic forces given by the gradients of T, p, and 
uP. Substituting the constitutive equations into the conserved 
currents i*' ^ 9 - * 1961 ). we have the first-order hydrody¬ 

namics in the Landau-Lifshitz frame. The terms in the other 
lines are the new terms appearing in the second-order hydro¬ 
dynamics. The second lines denote the relaxation terms given 
by the temporal and spatial derivatives of the dissipative vari¬ 
ables, which describe the relaxation processes of the dissipa¬ 
tive variables to the thermodynamic forces. The third, fourth, 
and fifth lines are composed of the products of the thermody¬ 
namic forces and dissipative variables, among which we re¬ 
mark that the vorticity term appears. The final lines give the 
non-linear terms of the dissipative variables. 

Our approach is based on a kind of statistical physics, and 
thus give microscopic expressions of the transport and relax¬ 
ation coefficeints. Here we present the resultant microscopic 
representations of the transport coefficients , i.e., the bulk vis¬ 
cosity C, thermal conductivity A, and shear viscosity 77 , and 
some of the relaxation times and 


C = -^(n,L-ifi) = C''°, 

\ = ^{jp,L-^jp) = \^^, 


m 


Tj 




(n, L ^n) _ oQ 
(n,L-in) “ " ’ 

{JP,L '^Jp) _ 
{Jp,L-^Jp) ~ ^ ’ 

{tP<^ ,L-^Tprj) 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


We leave the microscopic expressions of other coefficients in 
Appendix |C] We first note that A^^, and 77 ^*^ are per¬ 
fectly in agreement with those of the Chapman-Enskog (CE) 
expansion method iFt^I . which we denote as A^^, and 
77 ^^. Here it is noteworthy that our expressions of the trans¬ 
port coefficients can be nicely rewritten in the form of Green- 
Kubo formula ifTSl - f^ in the linear response theory. To see 
this, we first introduce the “time-evolved” vectors defined by 


(np(s),i;(s),^p^‘'(s)) = 



( 106 ) 
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where the time-evolution operator is given by the linearized 
collision operator. Then, we have 

1 ~ 

= ds(n(o),n(s)), (107) 

1 

= MJ^{0),Ms)), (108) 

1 

^ lOT io 

We note that the integrands in the formulae have the meanings 
of the relaxation functions or time correlation functions; 

i?n(s) = ^(fl(0),n(s)), ( 110 ) 

Rj{s) = -^{J^{0),Ms)), ( 111 ) 

^^(s) = Y^(d-''‘'(0),d-;.i.(s)). (112) 

We stress that the results of the transport coefficients all show 
the reliability of our approach based on the doublet scheme in 
the RG method. We remark that the naive version of moment 
method by Israel and Stewart (IS) fails to give the Chapman- 
Enskog formulae 12^ . as is well known. 

Thus it may be a good news for us that the explicit formu¬ 
lae of the relaxation times given above also differ from those 
given by IS ll^ . which read 


IS ^ (n,n) 

" (n,Ln)’ 

.J.IS _ { J^i) 

' ~ {jp,Ljpy 
^is ^ 

{'kP'^ ,LtTp„) 


(113) 

(114) 

(115) 


Indeed we shall now show that our formulae of the relaxation 
times allow a natural interpretation of them. To see this, we 
rewrite the expressions of the relaxation times given in Eqs. 
(I103b - (ll05b in terms of the time-evolved vectors again: 


,RG 


^RG 


^RG 


/“ds si?n(s) 
/“dsi?n(s) ’ 
dssRj{s) 
fo°^dsRj(s) ’ 

fy'dssR^(s) 
/(T dsi?^(s) ■ 


(116) 

(117) 

(118) 


It is noteworthy that all the relaxation times are expressed in 
terms of the relaxation functions i?n(s), Rj{s), and i?,r(s), 
respectively. Then the formulae (II 16b - (ll 18b allow the natural 
interpretation of the resepective relaxation times as the corre¬ 
lation times in the respective relaxation functions. We empha¬ 
size that it is for the first time that the relaxation times are ex¬ 
pressed in terms of the relaxation functions in the context of 
the derivation of the second-order relativistic hydrodynamic 
equation from the relativistic Boltzmann equation. 


C. Discussions 

We now examine the basic properties of the resultant hydro- 
dynamic equations ( l95l l and (l97])-(l99b. Eirst we show that our 
equation is really causal in the sense that the velocities of any 
fluctuation around the equilibrium is less than that of the light 
velocity with a detailed proof is left to Appendix |Dl where 
the stability of the static solution is also prooved. Next we 
compare our formulae of the relaxation equations with those 
derived by the moment method. Then we give numerical re¬ 
sults of the transport coefficients and relaxation times given by 
Eqs. (I100b - (ll02b and (I103b - (ll05b . respectively, and compare 
them with those by other methods. 


1. Causal property of hydrodynamic equations obtained by RG 
method 


We give a brief account of the proof that the velocities of 
hydrodynamic modes described by the hydrodynamic equa¬ 
tions ( |95] ) and (I97li-(l99b do not exceed the speed of light, i.e., 
the unity. We note that the detailed proof is presented in Ap- 
pendixlDl 

Eirst, we linearize the hydrodynamic equations around 
equilibrium state specified by constant temperature, constant 
chemical potential, and constant fluid flow, as follows: 

(A SX^siA ; A:) = 0, (119) 

where the matrices and are defined in 

Eqs. (ID23b - (ID26b and (ID33b - (ID36b . respectively, and the 
variables SXap (A ; k) are Fourier-Laplace transformations of 
SXa/ 3 {T; a) given by 


SXp, = 

SXpi = 
SXip = 


< 5 X 44 = 


A 




3TC^G 

h 


sn- 


eq 








eq 


yp2 ;^RG 

1 


< 54 , 


T 


eq 

5Ua 


eq 


J^2 


<5r, 




eq 

5T, 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 


with the arguments (r ; a) being omitted. Here, 5T, Sp, Su^, 
i5n, SJP, and are fluctuations from the equilibrium state 
and all coefficients take values at the equilibrium state. Fur¬ 
thermore, iA and k^ are frequency and wavelength conjugate 
to T and aP, respectively. We note that k^ is space-like vector, 
k^ < 0 , for any k^ 7 ^ 0 , which satisfies k^ = X^'^k^ because 
of aP = AP'^a^. We also note that the condition SX 7 ^ 0 into 
Eq. (II 19b leads to the dispersion relation A = A(fc). 

Then, as a typical quantity used for the check of the causal¬ 
ity, we examine a character velocity Uch that is defined as 


Uch = lim 

— —^■OO 




( 124 ) 
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With the use of the explicit definitions of and 

we can show that 

t^ch < 1, (125) 

is satisfied for any collision operator Lpg, that is, any differ¬ 
ential cross section. We emphasize that our hydrodynamic 
equations surely have the causal property, and hence can be 
applied to various high-energy hydrodynamic systems. 


2. Relation between relaxation equations by RG method and 
those by the other formalisms 


The relaxation equations (I97li-(l99]l can be made into the 
different form by iteration. Here, let us focus on the relax¬ 
ation equation for the stress tensor, i.e., Eq. (l99T l. by setting 

n = J'" = 0; 


d 


= 277^°TTp^ -f h 
dr 






(126) 


By solving this equation with respect to in an iterative 
manner and using the equality 


or or or 




dr 3 


Oa^'' 


(127) 


and the balance equation (l95T l. we find that the resultant equa¬ 
tion includes the following terms 

5 rjn f rj^ rj^ 5 

(128) 

in addition to those given in Eq. (I126E In this iterative man¬ 
ner, our hydrodynamic equation apparently gets to have all 
the terms given by of Eq. (73) in ll4^ . Notice, how¬ 
ever, that the last two terms of Eq. ( 1128b have a form of the 
second-order spatial derivatives of hydrodynamic variables, 
which make the hydrodynamic equation parabolic and accord¬ 
ingly acausal. Hence, we have an important observation that 
the naive iteration may spoil the causal property of the orig¬ 
inal hydrodynamic equation, and thus we must use the origi¬ 
nal form of the relaxation equations (1126b or (I97b-(l99b. Eur- 
thermore, since the appearance of the nonlinear vortex term 
^ seems to be inevitably associated with that of the 
second-order spatial derivative terms, the explicit appearance 
of such a nonlinear vortex term should be avoided in the re¬ 
laxation equation although its effect should be included in Eq. 
(1126b implicitly. 


3. Numerical example: transport coefficients and relaxation times 

In this subsection, we present numerical examples of the 
transport coefficients and relaxation times using the micro¬ 


scopic expressions given in the present approach, and com¬ 
pare them with those in the previous works. Note that the 
microscopic expressions are solely given in terms of the lin¬ 
earized collision operator L, which is in turn uniquely deter¬ 
mined by the transition probability a;(p, pi|p 2 , Pa)- A gen¬ 
eral form of the transition probability reads 

w(p, Pi|p 2 , Pa) = 5^(p-f Pi -P 2 -P 3 )sa{s, 9), (129) 

where a{s, 9) denotes a differential cross section, s = {p + 
Pi)^ a total momentum squared, and 9 = cos“^[(p — pi) ■ 
(p 2 — Ps)/{p — Pi)'^] a scattering angle. Here, we examine the 
case of a constant cross section for simplicity; 

cr(s, 9) = (130) 

with (Jt being a total cross section. 

We focus on the shear viscosity p and relaxation time 
for the stress tensor in the classical and massless limits, i.e., 
a = 0 and m/T — 0. The calculation of p and can be re¬ 
duced to that of Xp = which satisfies the linear 

equation The last equation can be solved nu¬ 

merically in an exact manner without recourse to any ansatz 
for the functional form of Xp. 

In Table |I] we show the numerical results together with 
those of the previous works. We confirm that our formulae 
for p and give results different from those by the (naive) 
Israel-Stewart moment method ll^ . Eurthermore, our relax¬ 
ation time differs from that of Denicol et al. which is 

an improvement of the Israel-Stewart moment method adopt¬ 
ing 41 moments, although their result for the shear viscos¬ 
ity tends to numerically agrees with the Chapman-Enskog/RG 
value 


TABLE I: Values of the shear viscosity and relaxation time for 
the stress tensor for a classical gas with a constant cross section in 
the massless limit, in the RG method, Israel-Stewart’s 14-moment 
method 12^ . and Denicol et al.’s 41-moment method 1^ . 


RG 

Israel-Stewart Denicol et al. 

p [T/ctt] 1.27 

1.2 

1.267 

Ttt [l/noT] 1.66 

1.8 

2 


IV. SUMMARY AND CONCLUDING REMARK S 


In this paper, we have derived the second-order hydrody¬ 
namic equation systematically from the relativistic Boltzmann 
equation with the quantum statistical effect. Our derivation is 
based on a novel development of the renormalization-group 
(RG) method. In this method, we have solved the Boltzmann 
equation faithfully in a way valid up to the mesoscopic scales 
of space and time, and then have reduced the solution to a 
simpler equation describing the mesoscopic dynamics of the 
Boltzmann equation. We have found that our theory nicely 
gives a compact expression of the deviation of the distribution 
function in terms of the linearized collision operator, which 
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is different from those used as an ansatz in the conventional 
fourteen-moment method. In fact, in contrast to the ansatz in 
the fourteen-moment method, our distribution function pro¬ 
duces the transport coefficients which have the same micro¬ 
scopic expressions as those derived in the Chapman-Enskog 
expansion method. Furthermore, new microscopic expres¬ 
sions of the relaxation times are obtained, which differ from 
those derived in any other formalisms such as the moment 
method. We have shown that the present expressions of the 
relaxation times can be nicely rewritten in terms of the re¬ 
spective relaxation functions, which allow a physically natu¬ 
ral interpretation of the relaxation times, and thus assert the 
plausibility of our results. 
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Appendix A: Detailed derivation of the excited modes and their 
explicit expressions 


The present asymptotic analysis utilizing a perturbation 
theory is based on the physical assumption that only the spa¬ 
tial inhomogeneity is the origin of the dissipation, and the ex¬ 
pansion parameter e is introduced for characterizing the inho¬ 
mogeneity, which may be identified with the Knudsen num¬ 
ber: The inhomogeneity gives a deviation of the distribution 
function from the local equilibrium one f^'^, and accordingly 
the ratio of the deviation to is necessarily proportional to 
e. We emphasize that the inhomogeneity and the ratio of the 
distribution functions are necessarily of the same order in our 
asymptotic analysis. It is worth emphasizing that the present 
asymptotic analysis combined with the perturbative expansion 
successfully solves the Boltzmann equation consistently and 
leads to the mesoscopic dynamics including the constitutive 
equations that relate the dissipative quantities and the spatial 
gradients of the equilibrium quantities. 


In this Appendix, we derive the expression of Eq. 
whose calculation can be reduced to that of 


[QoFo]p - [Fq - PqFo]^ - Fop - y^opVoapiy^OjFo), 

(Al) 


with 


Fop — 


p ■ u 




T 


(A2) 


Here, we have used Eq. (l44l i and F}/t _ q] 

We introduce the following quantities for later covenience: 


a, = j dprp^f;^ip • uY, f = 0, 1, • • • . (A3) 


We have given a proof that the propagating velocities of the 
fluctuations of the hydrodynamical variables do not exceed the 
light velocity, and hence our seconder-order equation ensures 
the desired causality. We have also proved that the equilibrium 
state is stable for any perturbation described by our equation. 
These results strongly suggest the validity of our formulation 
based on the RG method. 


Then the metric = {ipQ, ipY) are expressed as 


Vo'' = asu^u'' + {m ai - , 


u4 4u n 

% = Vo = 

44 

Vo =«!. 


(A4) 

(A5) 

(A6) 


We have demonstrated numerically that the relaxation times 
differ from those given in the moment methods in the litera¬ 
ture, even in the sophisticated one so as to numerically repro¬ 
duce the transport coefficients given in the Chapman-Enskog 
(and RG) method. The calculation was done only in the case 
of a constant differential cross section. It is interesting to ex¬ 
tend the present calculation to more realistic cases with the 
differential cross section depending on the total momentum 
and scattering angle, which have immediate applications to 
relativistic systems consisting of quarks, gluons, and hadrons. 
Then it is an imperative task to apply the present method to de¬ 
rive the multi-component relativistic hydrodynamic equation, 
which is now under way IsHl . It is of interest to use the re¬ 
sultant equations for phenomenological analysis of relativistic 
heavy-ion collisions performed in RHIC and LHC, although 
there exist multi-component hydrodynamic equations deriven 
on the basis of the moment method 182110 . which was found 
to have unsatisfactory aspects for the single-component equa¬ 
tion, as was shown in the present article. 


while the inverse metric Voap 


VOfiv ~ 


aiu^u'^ 


+ 




VotJ.4 — Voip 


ttaoi — 02 w?ai — 03 
—a2U^ 


0301 — 02 ’ 


7o44 — 


as 


03O1 — O2 


The inner products ((^g, Fg) are evaluated as follows: 


^ m'oi - 03 

Wo,Fo) =--- 


2^2 2^ 


_ TO" Op - 02 

3 F’ 

4 TO^Oo - 02 1 

{(Po,Fo) = - - --u. 


(A7) 

(A8) 

(A9) 


(AlO) 

(All) 
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Inserting the inverse metric (IA7l i-( IA91 l and the inner prod¬ 
ucts dAlOb and (lAl lb into Eq. (lAlb . we have 


[QoFo]p 


1 1 
T p - u 


np(-v • u) - J^T 


m?a[) — 02 
m?ai — 03 



with s being the entropy density. We note that the relations 
(IA25b and (IA26b can be shown not only by the Gibbs-Duhem 
equation but also by a straightforward manipulation based on 
the explicit forms of n, e, and P. 




(A12) 


Appendix B: Solution to the linear differential equation 
with a time dependent inhomogeneous term 


where lip, J^, and are given by 

m2(o20o-of) 2 

np = —57- Trip • 

3(0301 - 02 ) 


-f 


771^(0300 — 0201), , 771 * 


3(0301 - O 2 ) 


-{p-u) - 




3 

[P'u) - 2 - 

m^ao — a2 


777^01 




As is shown below, the following relatons hold; 

777^(0200 - of) 


3(0301 - of) 


-i dP dn _ dP dn 

£ _ dT d(i dfi dT 

o de dn _ de dn 

dT dfi dfi dT 


1 


dP de 


9 / \ c/c 

777^(0300-0201) _ dTdJl 


dP de 
dfj. dT 


3(0301 - of) 


dn de 
dfi dT 


dn de 
dT dfi 

2 dP 

Qi - Q3 _ ^ £ 

m^ao - 02 ^ 


dn 


(A13) 

(A14) 

(A15) 

dP 

(A16) 

(A17) 

(A18) 


Then we see that [Qo^ojp in Eq. (IA12b takes the form given 
in Eq. (l58l l. In the derivation of the above relations, we have 
used the equations derived from the explicit forms of n, e, and 
P given by Eqs. (I22l i- (l24l i. 


_ 1 
^ = “^“2 


dn 

Tr 

dn 1 


A 

rp2 


Oi, 


de 


1 

J'2 

1 


O3 


2^2 


02, 




dP 


1 


Tr ^ ^ ~ Apoo), 

dP 1 . 21 


(A19) 

(A20) 

(A21) 

(A22) 

(A23) 

(A24) 


In this Appendix, we present the detailed derivation of the 
second-order solution dSOl l and initial value dSTI) . We rewrite 
the second-order equation (iTSl l into 

^A(t) =LX(r)+iT(r-ro), (Bl) 

withA(r) = (rVA“'/("HT). 

The solution reads 

X{t) = e^^'^“’^“A(To) -f [ dT'e^^'^~^ '‘Kir' — tq) 




X{to)+ dT'PoK{T'- To) 


+ f '^'^QoK{t' - To), 

J Tn 


(B2) 


where we have inserted 1 = Pg + Qo in front of K{t' — tq). 
Substituting the Taylor expansion 


K{t' - To) = 


s=0 


(B3) 


into Eq. (IB2b and carrying out integration with respect to r', 
we have 


A(r) 


= qL(t-To) 


XiTo) + 


_ pp--ro)d/ds\j_ 


(1 -e' 


)i-d/ds)-^Po 


+ _ ^(r-ro)d/ds^^j^ _ 


= e 


L{t-to) _ p(r-ro)a/as 
Lir-To) 


K{s] 


s=0 


— e 


A(to) +Qi{L- d/ds)-^QoK{s) 

(1 - e(^-^'>)^/^")(-5/9s)-Ao 

L(t-to) _ p(r-ro)a/as 
{'r—To)d/ds tf —1 


s=OJ 


+ - e''^-^°^‘=>/^‘^)Pi{L - d/ds)-^Qo 


^Qi{L-d/ds)-^Qo 


K{s) 


s=0 


(B4) 


and the relations derived from the Gibbs-Duhem equation 
dP = sdP + ndp. 


dP _ _ e + P — pn 

T 
dP 


(A25) 

(A26) 


where 1 = Pq + Pi + Qi has been inserted in front of 
{L — d/ds)~^QoK{s) inthe second line of Eq. (IB4b . We note 
that the contributions from the inhomogeneous term K (r- tq) 
are decomposed into two parts, whose time dependencies are 
given by and ed'-'ro)d/ds^ respectively. The former 

gives a fast motion characterized by the eigenvalues of L act¬ 
ing on the Qo space, while the time dependence of the latter is 
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independent of the dynamics due to the absence of L. Since 
we are interested in the motion coming from the Pg and Pi 
spaces, we can eliminate the former associated with the Qi 
space with a choice of the initial value X (tq) that has not yet 
been specified as follows: 


^(ro) = -Qi{L - d/ds)-^QoKis) 

which reduces Eq. (IB4t to 


s=0 


(B5) 


X{t) = f(l - 

+ (e^(—o) _ - d/ds)-^Qo 


— Q{'r-'ro)d/ds 


Qi{L-d/ds)-^Qa K{s) 


s=0 


(B 6 ) 


Using X(to) = we can convert Eqs. (IB5b 

and (IB 6 I 1 into Eqs. dsTT i and dSOl l. respectively. 


Appendix C: Detailed derivation of the relaxation equations 

In this Appendix, we present a detailed derivation of the 
relaxation equation given by Eqs. (193-®. 

Eirst, we introduce the differential operator given by 


where 


1 1 d 

(P • ^ + ep • V 


^pq (P ■ '^)'^pn^o 


(Cl) 


D„ = 


= 

^pq — 

^ Pi'^pq'i ^ — Pi 

p • U 

(C2) 

Spq^ 

a = 4, 


eV^, 

a = p, 


d 


(C3) 


a = 4. 



Then Eq. (I93b is converted into the following form: 

+ [i-'xVfe]) +0(e3), (C4) 

where we have introduced the following “vectors” consisting 
of three components: 

= (C5) 

~ (n, t/p, TTpcr)? (C6) 

= (n, J^.TT^x), (C 7 ) 

Xi = iflp/{-TC^%hJP/iT^X^^),7TP//i-2Tr,^^)), 

(C8) 

(C9) 


with i, j, and k being indices specifying the vector compo¬ 
nents. 

We expand the left-hand sides of Eq. (IC4b as 


+e{L-^\ [nni-^x^] 

+ e{L-^\v^L-^X^)D^P^,. (CIO) 


The first and third terms of Eq. (1C 10b are calculated to be 


{L-^\ {nn)-\^D^n) = e{L-^\x^)X’, (Cll) 
e{L~^'ip\v°'L~^x’)Da'ipj 

(C12) 


respectively. In Eq. (ICl lb . we have introduced 


Xl = ■u,X^^T^h-^Vp{^^/T),2rJ^^XpU,). 

(C13) 


Substituting Eq. (1C 10b combined with Eqs. (ICl lb and (IC12b 
into Eq. ( IC4b . we have 


e{L ^■{p\x^)iPj 
= e{L-^\ x^)X' + e{L-^\ 

+ g 2 ^ + (9(e3), (C14) 


with 


M 


i,j,k _ 








(C15) 


[nnL-^x^\). 

(CI 6 ) 


Some coefficients in Eq. (1C 14b can be easily calculated as 





































17 


(L-V\X'') = 


/ 

(L-in,n) 


{L-^fl, ttP"^) 

\ 



_y^RG 

y2;i^RG 

-2TpRG 




(L-ij^,n) 


{L-^JP,ttP^) 




_y^RG 

y2;^RG 

-2TpRG 





h{L-HP'',jP) 

{L-^ttP'^ ,ttP'^) 



V 

_y^RG 

y2;^RG 

-2TpRG 

) 


/ 

(L-in,L-i 

n) h{L-^fi,L 

-i>) (L-i 

n, 

L- 


/ 1 0 0 

0 hAi^P 0 

^ 0 0 


{L-^\L-^X^) = 


-rc^° _27 ’^rg 

(L-1>,L-In) 

r^RG y2^RG 27’^RG 

h{L-'^7tP'',L-^JP) {L-^ttP'',L-^ttP'^) 


{L-^\v<^L-^X^) = 


_2^^rg 

2^2 ^RG 

-2TpRG / 

1 

o 

“ 'l 


0 -hrjAPP 

0 , 


1 

o 

o 

^/iivpcr / 


(L-in,u“L-in) 


(L-ilT,u“L-i7r^'^) 

_2^(^rg 

y2^RG 

-2Tr]^G 

{L-^Jp,v^L-^7rP'^) 

_y^RG 

2^2;^RG 

h{L-^nP’',v‘^L-^JP) 

_2TpRG 

_2^^rg 

y2^RG 

-2Tr]^G 


V 

/ 0 -WnjA“P 0 

-fjnA'"® 0 -tj^AP'^P^ 

0 -M^jAP^^P 0 


Here, we have introduced the relaxation times rn, tj, and 
and the relaxation lengths £uj, ^jn, and They are 
defined as follows, 


r 


^jn = 




_ {L-^Ii,L~^A) 

''’n — y^RG ’ 

(C20) 

- sT^A^g 

(C21) 

(L-iTf'^AL-iTf^,) 

“ lOTpRG 

(C22) 

^ _ {L-^U,vPL-^X) 

^nj — 32^2 ;^RG ’ 

(C23) 

1 

{An = l^^XAX), 

(^A^pAp<^^ — 

5 




3TC^G 

JP,v''L 
lOT^ 

{L-^7rP‘',v^L-^X) 


“ lOTTyRG 


(■ttJ = — 


5T2ARG 


(C17) 


(CIS) 


(C19) 


(C24) 

(C25) 

(C26) 


•RG Tj®, and 


We note that rn, tj, and are denoted as Tj 
in the text and given in Eqs. (I1031 l- (I105I) . 

From now on, we examine the terms associated with 
and N'’’^ in Eq. (IC141 l. We write down the useful formulae for 
space-like tensors A for later convenience; 


(C27) 

(C28) 

(C29) 


= lAP'^AP^iA^^Xp) + ^-AP^P^A^<^P\^p)) + 

<30 ' '' 


(C30) 
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_|_ ^^{fj,v) (per) {a/3) _|_ (po-> (a/3) ^ 

3 


15 

I- 


35 

«> 


(C31) 


where we have defined = i(A'^^A‘'‘^ — A^'^A'-'p) and 

AiP’') = VtP'^P'^Ap^ for an arbitrary tensor A^''. In the first 
equality of Eq. (IC291 l and the first equality of Eq. ( IC3U , we 
have used the fact that a space-like rank-two tensor 
is decomposed to be 

BP'' = AP''BP pji + B''P'''> + B^P''^. (C32) 

The numerical factors may be verified by contracting both 
sides of equations. To see how to use these formulae, let us 
consider {A^P'^'>°‘^)ip(p^)Xai 3 , which is found in the fifth line 
after the first equality of Eq. ( IC49l l. for instance. 

= (C33) 

where we have used Eq. (IC28l l in the second equality. 

Using the formulae (IC271 l. (IC28b . and (IC30b . the non-linear 
terms of for = H, JP, ttP" can be reduced to 

=fl 




= —e 


— e 


(L-in,s[L-in,L-in]) 

2{Tcy 


2{T^X/hy 




JpJ(T 


2 (L-^n, L-iTf"'^]) 


= e-\ - 


2{2TriY 




J'^Jo 


V 2(T02 

6(r2A/fi)2 
io(2rp)2 

= e^(fennnn^ + bujjJPJp + bmnr'^^'^T^ptr), 


(C34) 


= JP 




= e 


.{L-^Jp,B[L-^U,L-^JP]) 

(rc)(T2A/fi) 


nj. 


^{L-^Jp,B[L-^JP,L-H°'/^]) 
(T2A//i)(2Tp) 
'{L-^JP,B[L~^fl,L-^Jp]) 


UJP 


JPi 


' 3(rc)(r2A/fi) 

b{T^\/h){2Tr/) 

= ^^{bjiijAJP + bjj-^JPjTpP), 
3/>* = 


-^-{L-^nP-',B[L-^^,L-^'^])xm 


, {L-^ttP'', B[L-^fl, L-^ttP'^]} 


flTT. 


— e 


(rc)(2Tp) 

2 {L-^nP'',B[L-^JP, L-i>]) 


2(r2A/fi)2 

L-i7r“^]) 

2(2rry)2 


JpJa 


pa '^OL^ 


(C35) 


5(TC)(2Tp) 

10(T2A//i)2 

(35/6)(2T77)2 p) 

= e^ib^n-.A^P'' + b^jjJ^Pr') + b^^^irP^Pir''^p), (C36) 

where the coefficients 6nnn, bnjj, bn^„, bjnj, bjj„, b^n^r, 
b^rjj, and b^„^ are given by 


_ B[L-^n, L-^fi\) 

fonnn = 2{TC)^ ’ 

_ {L-^fl,B[L-^JP,L-^Jp]) 

bnjj = - 


bllTTTr — 

bjnj = 

bjjTT = 


Q{T^X/hf 

10(2rr7)2 
_ {L-\jP,B[L-^%L-\jp\) 
3(TC)(T2A/fi) 

5{T^X/h){2Tr]) 


(C37) 

(C38) 

(C39) 

(C40) 

(C41) 
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“ 5(rC)(2T?7) 

^JJ - 10(T2A//i)2 

(L-l7^^^i?[L-l7^;i,L-l7^A.]) 

(35/6)(2T7?)2 

Next, we rewrite TV® -' as follows: 


(C42) 

(C43) 


(C45) 


(C44) 


The temporal derivative of T, fi/T, and are rewritten by 
using the balance equations up to the first order with respect 
to e, which correspond to the Euler equation: 




— +evV 

OT 


rrL-^xni’j 




d 


+ {nr)-^vf^-^[rnL-^x^])^,vpT 


At=-T — 
dr de 


e9 + 0{e^), 


d /r 


dP 

dn 


ee + 0{€^), 



leV'®T+^eV'®^ + 0(e2). 


(C46) 

(C47) 

(C48) 


Using the formulae (IC27l) - (IC311 l and Euler equation (IC46b - 
( IC48I ), we convert Eq. ( IC45I ) into the following forms: 


= n 


— + evV 

OT 


dP 




d dP 

'Wf~ 


e '^T 


d Id 




e(^L-% {nnr^ 

+e"(L-in, (nnr^ 

J J _2Tn / p 


dT Tduii\ T^X/h 


d T d 


/eq/eUL-ln 

~c 

JpVpT 


m 


5f hdup\ T^X/h 


feqfeqp-ljp. ^ 




= e 


/ ^ ^ p 

/L-^n, -T 

A“/5 / 

^_^^-ln^(/eqJeq)-l 


-l^po- 

-2T^ 
d dP 

'Wf~ 


de 

d 1 d 


e ^ T -i 


/eq/eq^-^n 


-TC 


m 


dT Tdu°‘ 
d T d 


h du° 


feqjeqp-lj^ 

T^X/h 

feqjFeqL-lj^ 


T^X/h 


JPVpT 


J^vp- 

H rp 




^afi-yS 


L-'n, innr^u 


, d nrpL-H^s\ 


du9 -2Tr] / 


}ttP'^(t, 


= e 


Knnm + J^V^T + K^uU^^P^ + ^n^Tr'^'^a 


JJJ^ ^ P rji 


(C49) 
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= > 


e(L-v>, innr^ i + nnL-^r)xj 


+ " \ ^ > [^ dT^Tdu^l -TC / 

+ " / J aT + ra^fl -2T77 


nv«T 


TTpcrVpT 


\ ’V. . / [ ar raw^J -2 t?7 / 

j/,-i™ ,,.<,r.,^-iL» a r a 1 nnL-‘ir^y_ „ 




- p2 t^/T-^T I'feqfeq'l-l _T — —-— — 


ap| a api a 1 


[ 3 \ p’w a y 1^ ae|„ar an|^afj t 2 a//i 

[,4 4 ^]/!!£^)nv.. 

- (/•’/-)- [..4 .14] 

[. 444 ] /44n)nv-| 


Aa/37(5 . 

^ /r-1 


-if a T a 1 nnL-^ir^s 


(/-/-)- + -^J 


A^'^A“^ / 






^afSjS 


L-^j^,inn)-^v. 


' “ au/^ r2;^4 

a nnL-^j0 


5 4 ^ ^ T2A4 r p 

^ , a ' 

+ 3 4 j ^75^5 ^24;^ 4p 


1p^ = TT^ 


e(L-'ir"7(/''i/“i)-‘ J]: + e»'V Hf'^L-'ipXi 


14 - SI. A] 


. .(.-.*.7 [.4 441 A44)a,v.r 

+ e^(L-4'^^ (/eq/eq)-i4^/4g^^nv^u 


r2A//i /'^p^Pt 
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^paOL 0 


du^ -2Tr] / 

,-eq,_ir ^api d dp\ d]nnL-^7:^p 


A pcra/3 . 


T-i- (\_T — — _— __ 

^ ^ [ de ^dT dn ^d§\ -2Tr^ 


' duf< -TC 


5 r ^ ^ r9r + T9u“J T2A//i r ^ ^ 

+ 7r,.,(/ / ) ^TT 

+ (/eq/eq)-i^^^ )7rP<^a->p 

= e" . (C51) 

The coefficients Knn, 'tn j, j, ttnir, ttj^, Kj^, K^pj, Kjj, Kjj, Kj^, Kj^, Kt^h, and n'^l are defined by 

™“\^ >(/ / ) [ ^ae„9r an^5f“^3 -TC /’ 

«(i) = fr<tf-en)-i L A + iA1 nnL-\i \ 

nj- 3 A ^ > ^dT^Tdui^ T-^X/h /’ 


^\_P<P, ,1^) 


ae|„9r an|^5f 3 au'^J -TC 

g 1 a 1 jeqjeq^^-lj^. 


-lU- 3 j.2X/h 

4..-(p-.n,,rv-.-[A-?^]AA 


a 

1 

a ■ 

yeq jeqj^- 

-iH' 

dT 

+ y 

du^^ _ 

-TC 


a 

T 

a ■ 

yeq Jeq^- 

-A 

af 

^ h 

au^_ 

-TC 


'T' 

dP 

a 

api 

a 

— 1 

de 


' an g a^ 


^m = ^U-'j innrA-T?^ 2._2£ 
^ AAa-ij ^ nnL-^l 

5 \ ^ ^ T2A/fi 




L-^jp, inn) 


_1 g jeqJeqj^-lj^ 

T2A/fi 


(rqfeq)-iL g , 1 ^ 1 

Jtt 5 \^ mA/ / 1 AaT^TSM^'J -2rr/ 

,) ^ aaA- iJ frvA-1 La, Z±.] nnL-n^ 

5 A f^^ij J > hdu''\ -2T7J 


_AA.^P-/^_1^ ,PeqJeq^-l,. a nni-^n 
Knn- g (T 7r^.,(/ / ) 



22 


,,, , 


7 = 

«S-^(i-v. (/•’/•’)- 

\^ivpa . 

/ 7—1 


9 1 9 

A ^ ^ 

^ h^p 


^(1) = 

- 5 




-T 


dP 

de 


nnL-^j„ 

T^\/h 

nnL-^ja 

T'^Xjh 

d 


1 ^ 9 


■(J2 = 


d dP 
dT dn 

d nnL-^T^pa 
duP -2Tf] 
d nnL-^^p. 


nni-^TTp^ 

-2Tt] 


duP -2Tr] 


(C64) 

(C65) 

(C66) 

(C67) 

(C68) 


Substituting the above equations into Eqs. ( IC141 l and set¬ 
ting e equal to 1, we arrive at the explicit form of the relaxation 
equations (I97b-(l99b. 


Instead of ST, Sji, and Su^ which are not independent of 
each other because 5u^ Uq = 0, we use the following vari¬ 
ables as the five independent variables composed of ST, h/i, 
and SuP', 


Appendix D: Proofs of stability and causality 


SXi^ = -Siu^/T) = -Su^/To +ST uop/T^,(DU) 
SXii = Sit^/T) = S^l/To - STf^o/Tl (D14) 


a. Proof of stability around the static solution 

In this Appendix, we show that the static solution, espe¬ 
cially the equilibrium solution, of the second-order relativis¬ 
tic hydrodynamic equation given by the pair of Eqs. ( |92] | and 
(l93t is stable against a small perturbation. 

A generic constant solution reads 


T(a; 

t) 

= To, 

(Dl) 

/r(cr; 

t) 

= Po, 

(D2) 

M^(cr; 

t) 

= 

(D3) 

n(cr; 

t) 

= 0, 

(D4) 

JP{(j-, 

t) 

= 0, 

(D5) 

TTP’^ia; 

't) 

= 0, 

(D6) 


where Tq, hq, and Uq are constant. We remark that the equi¬ 
librium state is correspondent to the special case of Uq = 
( 1 , 0 , 0 , 0 ). 

To show the stability of the constant solution, we apply the 
linear stability analysis to the second-order relativistic hydro- 
dynamic equation (l92l l and ( l93t . We expand T, /r, u^. If, Jp, 
and TtP'^ around the constant solution as follows: 


T(ct; r) 

= ro + hT(a; r), 

(D7) 

H{a; t) 

= Po + ^/4(cr; t). 

(D8) 

uP{a ; r) 

= Mq -1- 5uP{(j ; r), 

(D9) 

n(CT; r) 

II 

(DIO) 

JP{a-, t) 

II 

(Dll) 

TtP^{a ; r) 

= SitP^a^T). 

(D12) 


In the following, we suppress the subscript “0” in Tq, /io, and 
Uq. Eurthermore, we introduce the following variables 


sx^, 

5X^i 


-Afj_y hn/3 _ Sitpu 

(n, L-in) L-i V)’ 


hSJ^ 


(D15) 

(D16) 


which are expressed in terms of <511, SJP, and SttP'^. We 
treat SX^p = {SX^^, SXfj_ 4 , SX 4 ^, SX 44 ) as the fundamen¬ 
tal variables. 

Substituting Eqs. (ID7b - (ID101 l into the second-order rela¬ 
tivistic hydrodynamic equation ( |92] | and ( |9^ , we obtain the 
linearized equation governing SXap as 


, ^^q)^SX4p + , L-^ipf)^SX,p 

+ {v^,vPip^Q)XpSX4p 
+ {^^,vPL-^^f)VpSX,p 
= 0, (D17) 

(L-vr > 

+ (i-Vr, 

,vPL-\f)VpSX,p 

= , LL-^^f)SX,p. (D18) 


In the derivation of Eqs. (ID 17b and (ID 18b . we have used the 
fact that 


We assume that the higher term than second order in terms of 
ST, Sfi, SuP, (5n, SJP, and SitP'' can be neglected since these 
quantities are small. 


(D19) 

(5(^p) = [L-Vr]/^pa, (D20) 
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with 


aoi _ 

Vip = 


-A'"'" n. 


p ^ '^p j ^ 

a = 4. 


We can reduce Eqs. (ID 17b and (ID 18b to 

^a/3.75 ^SX^s + <5X^5 = 0, 

or 


(D21) 


(D22) 


where and are defined as 


j^pP,vS ^ (L-Vi'^, ), (D23) 

= (L-Vf ,‘^o), (D24) 

(D25) 

^4/5.4^ = (D26) 

gp,p,u5 ^ 

+ (L-Vf (D27) 

j^pPAS = (L-Vf (D28) 

= ((^^,u'’L-Vf)Vp, (D29) 

54/3.45 ^ ^^P^yP^l)Xp. (D30) 


We convert Eq. (ID22b into the algebraic equation, using 
the Eourier and Laplace transformations with respect to the 
spatial variable and the temporal variable r, respectively. 
By substituting 

5X^p{a ■ t) = 6Xcp{k ; A) e*''— (D31) 

into Eq. (ID22b . we have 

(A - 5“/3.^^) SX^s = 0, (D32) 

where is defined as 

+ {L-^<pf ,vPL-^^f)tkp, (D33) 
5/^/3.4^ = ,^p^g)ifcp, (D34) 

54/3.-^ = {p^^,vPL-^^f)ikp, (D35) 

54/3.45 ^ (D36) 


We note that k^ is a space-like vector satisfying k^ = A^'' k,y. 
In the rest of this section, we use the matrix representation 
when no misunderstanding is expected. 

Since we are interested in a solution other than 5X = 0, we 
can impose 

det(AA-5)=0. (D37) 

It is noted that Eq. (ID37b leads to the dispersion relation 


A = A(fc). 


(D38) 


close to the zero with time evolution. Therefore, our task is to 
show that the real part of A(fc) is positive for any k^. 

We show that A is a real symmetric positive-definite matrix 
as follows: 


WaP A“/^’'/''3 W^S 

= ( Wp,p 5“4 ppP + ipP ^ ^-1 (^(j'5 _|_ 

= (x.x)>0, Wap^O, (D39) 

with Xp = Wya [A“4 (/35‘“]p + W4a ^Sp- I" Eq. ( |D39b , we 
have used the positive-definite property of the inner product 

ii). 

Equation ( ID39b means that the inverse matrix A“4 exists, 
and A“4 is also a real symmetric positive-definite matrix. 
Thus, with the use of the Cholesky decomposition, we can 
represent A“4 as 


A-4=*5 5, (D40) 

where U denotes a real upper triangular matrix and ‘5 which 
is a transposed matrix of 5. Substituting Eq. (ID40b into Eq. 
(ID^ . we have 


det(A/-55‘5) = 0, (D41) 

where / denotes the unit matrix. It is noted that A(k) is an 
eigenvalue of 5 B*U. 

We find that the real part of A{k) is positive for any 
fc/' when Re(5 5*5) is a positive definite matrix where 
Re(M) = (M -f Mt)/2. In fact, we can show that 
Re(5 5 *5) is positive definite as follows: 

Wc,p [ReiU B*U)r^’'^^ w^s 
= Wc,p[URe{B)*U]°^l^’'^^w^s 
= [wUU[Re{B)r^^'^^[wUU 
= -[wU]pp , L-^ ^f)[wUU 

= -(V^, 5-4?/') >0, Wo^p^O, (D42) 

with ipp = [wU]pa ‘Pip - The inequality in the final line is sat¬ 
isfied because the vector ipp belongs to the Qo space spanned 
by the eigenvectors correspondent to the negative eigenval¬ 
ues of Lpg. Therefore, we conclude that the constant solution 
given by Eqs. (IDlb - (ID6b is stable against a small perturbation 
around the general constant solution. 


b. Proof of causality 

Here, we show that the propagation speed of the fluctuation 
5Xap is not beyond the unity, i.e., the speed of light. Here, 
we suppose that the propagation speed of SX^p is given by a 
character speed, whose Lorentz-invariant form may be given 
by 


Uch = ■\/-A 


P-u ‘"ch ‘'ch- 


The stability of the constant solution given by Eqs. dPTTi-dDfiTl 
against a small perturbation is equivalent to that 5X becomes 


(D43) 
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Here, we have introduced the space-like vector defined in 


terms of A(fc) given in ( ID38b as 


= lim 

— k'^—¥oo 




(D44) 


By differentiating Eq. (ID41b with respect to i we find 
that is an eigenvalue of U *U, i.e., 


det 




= 0 , 


with 


lim 


d r. 


CK^,'y 6 


(D45) 


(D46) 


— k^—¥oo 

whose components are given by 

^ , yP L-^ (D47) 


[CP 

[CP]^^’^^ = Wlvp^t)- 


(D48) 

(D49) 

(D50) 


An expectation value of U *U with respect to an arbitrary 
vector w' = *{U~^)w can be written as 




J Oc 0 


7(5 


r_1 \1 Oi' (3' ,'y'S' 


Wa'P'[U 7]' 'PPl’S' 

{x,vPx) 


tu, 


Wa'/3' A“'/5'.7'<5' wys' ( X , X) 


= {vP)^, (D51) 


with Xp = Wpa [L ^ -f W 4 a if Op- Here, we have intro¬ 


duced 


ioK^ 


{x, Ox) 
(x, x) 


(D52) 


with O being an arbitrary operator. 

It is important to note that if the inequality 


{VP)^{V'')^ < 1, 

are satisfied for any Xp, we can conclude 


(D53) 


= <1. (D54) 


Indeed, we can show that the inequality (ID53b is satisfied in 
this case. The proof is given as follows: First, with the use of 
the identities 


A yPy’' - 

P'' P P {p ■ U)2 

— < 1, (D55) 

(1), = 1, 

(D56) 

we obtain 


{-Ap,vPvn^<l. 

(D57) 

Then, we notice 


i-Ap^vPvn^ 


= -Ap, + { -Ap, SvP Sv'' 

>-Ap^{vP)^{vny 

(D58) 

where Syp^ = Syp Spq with SyP = yP — 

{yP)y., because 

- Ap^ Syp Syp > 0, 

(D59) 


due to the fact that Sv!^ is also a space-like vector. By combing 
Eq. (ID58b with Eq. (ID57b . we complete the proof. 

Thus, our fourteen-moment equation given by Eqs. (l92l i 
and (l9^ respects the causality in the linear analysis around 
the homogeneous steady state (IdT1i-(ID61i. in addition to the 
stability around the static solution. 
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